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^vq In this paper we wish to study hadronic wave functions using an instanton model for 



the QCD vacuum. The wave functions are defined in terms of gauge invariant Bethe 
Salpeter amplitudes which we have determined numerically using a Monte Carlo simula- 
tion of the instanton ensemble. We find that the pion and the proton, as well as the rho 
meson and the delta have very similar wavefunctions but observe a sizeable splitting be 

>: 

tween mesons or baryons with different spin. We compare our results with data obtained 
in lattice gauge simulations. 
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1. Introduction 

Hadronic correlation functions are an important tool in understanding the complicated 
nonperturbative structure of the QCD vacuum |T[. These correlation functions are defined 
as the vacuum expectation value Ii(x) =< T(jh(x)jh(0)) > of the time ordered product 
of two currents with the quantum numbers of a given hadron h. At short distances, 
asymptotic freedom predicts that the correlation function is dominated by the propagation 
of free quarks. In this case, the correlator follows a simple power law behavior, U(x) ~ 
1/x 6 for mesons and U(x) ~ 1/x 9 for baryons. Corrections to this behavior can be studied 
using the Operator Product Expansion. 

In general the correlation function can be expressed in terms of the contributions 
of intermediate physical states. At large spacelike separations, the lowest resonance in 
a given channel dominates and the corresponding mass and coupling constant can be 
extracted from the exponential falloff of the correlator. Roughly speaking, the coupling 
constant measures the probability for all quarks inside the hadron to be at the same 
point. In order to obtain more detailed information about the structure of the hadron 
the concept of hadronic wavefunctions was introduced in ||. In the case of the pion this 
quantity is defined as the following gauge invariant Bethe-Salpeter amplitude 

^(y) = [d 4 x <0\d(x)Pe l S* +VA{x ' )dx ' l5 u(x + y)\7i> . (1) 



Here the two quarks are connected by a string of flux. One should note that in a relativistic 
theory the concept of a hadronic wavefunction is not uniquely defined. Other functions 
providing somewhat different information include the Bethe-Salpeter amplitudes in a fixed 
gauge or density-density correlation functions (inside a given hadron). In practice the 
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Bethe-Salpeter amplitude eq. ([I]) is extracted from the three point correlator 

n w (x,y) = <0\T(d(x)Pe t C yA(x ' )dx ' l5 u(x + y)d(0) l5 u(0))\0>^i,(y)e- m ^ (2) 

where x has to be a large spacelike separation in order to ensure that the correlation 
function is dominated by the ground state and y is the separation of the two quarks in 
the transverse direction. 

The Bethe-Salpeter amplitude eq. ([D has been measured in a number of lattice gauge 
simulations, both at zero and at finite temperature 0, |[ £|, ||. However, with the exception 
of one work in the context of the bag model || and the case of very large temperatures, 
where the wave function in the dimensionally reduced theory can be evaluated from an 
effective Schrodinger equation 0, |8j , there have been no efforts to determine these wave- 
functions in a theoretical model. 

We would like to fill this gap by applying the Random Instanton Liquid Model 0, 



10| , pTTJI to the problem of calculating hadronic wavefunctions. The model is based on 
the assumption that the QCD vacuum is dominated by a particular type of gauge field 
configurations, small size instantons. The vacuum is then characterized by the density 
n = 1 fm 4 and the average size p = 1/3 fm of the instantons. Essentially, these numbers 
were chosen to reproduce two global properties of the QCD vacuum, the value of the 
quark and gluon condensates. Recently, we have determined a large number of hadronic 



two point correlation functions using this model [11, 12 1. The model not only successfully 
describes the masses and coupling constants of the lowest resonances in most channels, 
but also reproduces the results of a recent lattice calculation of two point correlation 
functions [ |T3]j . 

In the present work we want to extend this comparison to hadronic Bethe-Salpeter 
amplitudes. First of all, this will enable us to check whether the correlators calculated 



in [11], [12[ really correspond to the propagation of hadronic bound states. In addition 
to that, one can study whether the size and shape of the corresponding hadronic wave 
functions is similar to what was found in lattice simulations. In general, these are very 
non trivial questions. The Instanton Model takes into account only one very specific 
non perturbative effect among the many that determine the structure of QCD. Whether 
instanton induced effects really dominate the physics of a particular observable has to be 
checked in each individual case. 

Another interesting question concerns the differences between the wavefunctions of 
hadronic states with different quantum numbers. In a potential model approach, these 
differences are due to the gluomagnetic forces that give rise to a spin-spin interaction 
between quarks. These forces are absent in our calculations. However, there also exists 
a spin-dependent instanton induced interaction between quarks, which qualitatively de- 
scribes the mass splitting between different states. The question is whether it can also 
account for the differences in the size of the corresponding bound states. 

Finally we would like to address a few questions related to the interpretation of lattice 
results on hadronic Bethe-Salpeter amplitudes. In particular we want to study the ques- 
tions whether the observation of localized Bethe Salpeter amplitude uniquely reflects on 
the existence of a hadronic bound state and whether measurements of correlation func- 
tions at rather short distances on the order of one fermi (or even less) are sufficient to 
deduce properties of the ground state in a given channel. 



2. Wave Functions in the Random Instanton Model 



In the following we will study Bethe-Salpeter amplitudes for the pion, the rho meson, 



the nucleon and the delta resonance. For this purpose we have considered correlation 
functions involving the following currents 

f( x ) = d(x)^u{x), (3) 

j£(x) = d{x) lfi u{x), (4) 

jP( x ) = e abc (u a (x)C l5 d b (x)X a (x), (5) 

3aA x ) = e abc (2(u^x)C^d\x))u c a (x) + (u a (x)C^u b (x))d c a (x)). (6) 

Here q^ denotes a u or <i-quark spinor with color index a and spinor index a and C is the 
charge conjugation matrix. The nucleon current given above is a linear combination of 



the standard Ioffe currents [12[ . Our choice here is motivated by the fact that this current 



was used in the lattice simulations reported in M. We have done calculations using all 



six nucleon correlation functions introduced in |12j and will comment on the importance 
of the choice of current below. 

We consider the case of exact isospin symmetry so we have simply chosen the flavor 
content of the mesonic currents to minimize the number of contractions contributing to 
the correlation function. In the case of the nucleon we have used a proton current and 
defined the corresponding wave function by taking the two u quarks at the same point 
and measuring the correlator as a function of the separation of the d quark. Similarly, for 
the delta resonance we have taken the A + state and again measured the Bethe Salpeter 
amplitude as a function of the distance between the two u quarks and the d quark. 

The Bethe-Salpeter amplitude defined in eq. ([I]) involves a diagonal correlation func- 
tion of two pseudoscalar currents. The same information, however, can also be obtained 
from non diagonal correlation functions. In the case of the pion, we use the pseudoscalar- 



axialvector correlation function 

H,- ai Ax,y) = <0\T(d(x)Pe i i: +VA{xr)dx ' 1 M^ + y)d(0) 1 , l5 u(0))\0>, (7) 

while in the case of the rho meson, we have also determined the Bethe-Salpeter amplitude 
from the off diagonal vector-tensor correlation function 

U p _ T ^x,y) = <0\T(d(x)Pe t f: +VA ^' )dx ' 1 Mx + y)d(OKMmU>- (8) 

Similarly, for the nucleon one can use non diagonal correlation functions between the two 



different Ioffe currents [12fl . The advantage of using non diagonal correlators lies in the 
fact that these correlation functions receive no contribution from free quark propagation. 
For this reason the dominance of the ground state over excited states at large but finite 
separation is expected to be somewhat better for non diagonal correlators. In any case, 
comparing the Bethe Salpeter amplitudes determined from diagonal and non diagonal 
correlators provides a nice check on whether the measured wavefunctions really do reflect 
on the groundstate in the given channel. 

In addition to the channels considered here it would also be interesting to study the 
wave functions of other mesons. Of particular interest in comparison with the rho meson 
is the Oi meson. This channel, however, is difficult to analyze since the a\ correlation 
function at large distances is dominated by the pion |]J and the off diagonal ai-tensor 
correlation function turns out to be very small. Furthermore, it would be nice to demon- 
strate that one does not observe localized Bethe Salpeter amplitudes in channels in which 
(at least inside the model) no bound states are expected to form. In the Instanton Model, 
one such channel would be the scalar isovector (5-meson) channel. Again, the analysis 
turns out to be very difficult since the fact that the interaction is repulsive implies that the 
measured correlation function at large distances is very small so that the wave function 
cannot be reliably determined. 



The instanton ensemble in QCD can be generally described by a partition function of 
the type 

Z = / n i [dn i exp(-5' i )] exp(-S int ) H f=1)Nf deb(iD + imf) (9) 

where [dfli] denotes the integration over the collective coordinates (orientation, size, posi- 
tion) associated with the instantons, Si = 8ir 2 /g 2 (p) is the action of an individual instan- 
ton and Si n t is the instanton-instanton interaction. Since the evaluation of this partition 
function still constitutes a formidable problem we consider the even simpler Random In- 
stanton Model. In this model we assume that except for the size, which we keep fixed, the 
distribution of the collective coordinates is completely random. In particular, adopting 
this assumption we neglect the effects of the fermion determinant det(il) + imf), which 
is similar to "quenched" lattice simulations we are going to compare with. 

The Bethe-Salpeter amplitudes are then determined from the quark propagator in a 
given configuration. The pion correlation function eq. (0), for example, is given by 

I^foy) = <TiS ab (x + y)(Pe i C VA ^ dx ') bc S ca (-x)> (10) 

where the trace is performed with respect to the Dirac indices. The calculation of the 



propagator is explained in detail in [11]. Here we only outline the general strategy. The 
quark propagator in the field of a single instanton can be written as the sum of a zero 
mode contribution and the propagator due to non zero modes. Both pieces are known 
analytically. In the many instanton configuration we treat the zero modes exactly by 
numerically inverting the Dirac operator in the space spanned by the zero modes of the 
individual instantons. The non zero modes, on the other hand, are taken into account 
only approximately, by including the first term in a multiple scattering expansion. The 



resulting quark propagator is given by 

S(x,y) = J2 < x \<Pi >< <t>i\~ — f — r— I0J >< <t>j\v > +Snzm(xiV), (11) 

where <f>i(x) is the zero mode associated with the instanton / and 

Snzm(x, y) = S (x, y) + J2( S nzm( x i v) ~ s o(%, y)) (12) 

i 

is the propagator due to non zero modes. Here, Sq denotes the free propagator and S^ ZM 

is the non zero mode propagator in the field of an individual instanton. The non zero 



mode propagator also has to be corrected for the effects of current quark masses |l T[ • In 
this paper we take the light quark mass to be m u ^ = 10 MeV and the corresponding 
effects are small. The Schwinger factor Pexp(i / A^dx^) is calculated numerically for any 
given gauge field configuration by expanding the pathordered exponential as an infinite 
product 

Pexp(i / A^dx^) = n»(l + iA^x^dxi^) . (13) 

In practice this expression is approximated by a finite product where the stepsize dxi is 
determined by the magnitude of the local gauge potential. 

In general, the inclusion of the Schwinger factor eq. fll3|) is expected to give an impor- 
tant contribution to the measured wave functions, since it corresponds to an additional 
string type potential. The instanton liquid, however, does not produce confinement and 
no string potential is expected to appear. At large separations, the Schwinger factor 
decreases exponentially Pexp(i/ r Adx) ~ exp(— mr) with m being the mass renormal- 
ization of a heavy quark in the instanton background [O]. Using the 'sum ansatz' for 
the gauge field of the multi-instanton configuration we find m = 65 MeV, in agreement 



with estimates made in |14|| . Employing the more realistic ratio ansatz an even smaller 



screening mass m = 20 MeV is found. In any case, the insertion of the path ordered 
exponential will only lead to minor effects on the wave functions we are going to discuss. 

3. Results and comparison with Lattice Gauge Simulations 

Before we come to the results obtained, let us discuss the Bethe Salpeter wave functions 
obtained by using simplified models for the propagator. Since in all numerical simulations 
to date the separation between the sources is rather modest (at most on the order of 1 
fm) it is important to understand what the size of non asymptotic effects in the measured 
wave functions are. 

In the simplest case we consider a free fermion propagator 

and take the path ordered exponential to be unity. Here all gamma matrices and coor- 
dinates are taken to be euclidean. Inserting this propagator in the pion correlator eq.(|2]) 
we get[] 



Il*(x,y) = — - - ■ J sa ^-^.(1- -- + ...)■ (15) 



N c x 2 - y - N c 5y 2 

At distances r = Vr on the order of one fermi, this corresponds to a 'wave function' 



with a mean square radius^] J(r 2 ) = 0.89 fm, quite different from the asymptotic result. 
This effect becomes even more pronounced if one considers the propagation of a massive 
constituent quark 

S(x) = i D'(m q ,r) + m q D(m q ,r) (16) 



^ince we are discussing short range effects, we are using a symmetric configuration in which the 
quarks are located at transverse separations — y/2, +y/2. For asymptotic wave functions, this difference 
should not matter. 

2 The mean square radius is defined by (r 2 ) = J d 3 rr 2 \ip(r)\ 2 where the Bethe Salpeter wave function 
is normalized to one. 



where D(m, r) = m/ (4iT 2 T)Ki(mT) is the coordinate space propagator of a massive scalar 
particle. Taking the quark mass to be 300 MeV, we find a wave function with a size of 



J (r 2 ) = 0.79 fm at a distance of r = 1.0 fm. Even at a distance of 1.5 fm, the wave 



function is still fairly well localized with a radius of J (r 2 ) = 1.1 fm. Note however, that 
so far the appearance of a localized wave function is a purely geometric effect and has 
nothing to do with the existence of bound states. 

One should note that the fact that the Bethe Salpeter wavefunctions have not reached 
their asymptotic shape can not be detected from an analysis of the spacelike screening 
mass. The screening masses in this simple model are m scr = 2m q for mesons and m scr = 
3m q for baryons. There are only small power corrections to this behavior that are due 
the fact that the quarks can have some transverse momentum. We would also like to 
emphasize that considering the effects of finite quark masses on the Bethe Salpeter wave 
functions is very important at finite temperature where quarks develop thermal screening 
masses m scr = ttT. Around T c this mass is on the order of 500 MeV . If the wave functions 
are determined at distances x ~ 1/T, which is typical for current lattice simulations, 
masses will have an important effect on the measured Bethe Salpeter wave functions. 

After these remarks concerning non asymptotic effects we would like to study how 
instantons lead to the formation of bound states that manifest themselves in localized 
asymptotic Bethe Salpeter amplitudes. A simple one-instanton approximation that has 
been used [15| to analyze the qualitative effects of instantons on the short- distance be- 



haviour of correlation functions is based on the following expression for the propagator in 
a dilute instanton gas 



,j . 0o, z O)</>o,*(O) 



S(x) = S (x)+n / rf 4 z '^ : ru ' zv ' . (17) 



—%mr 



Here n is the density of instantons, (t) ,z{^) is a zero mode localized at z and m* 



2vrp c wno/6 is the effective quark mass in the instanton gas. Recall, that such propagator 
eq. (|17D leads to an attractive pion correlation function and explains the failure of the 
operator product expansion in this channel. At distances greater than r = 0.5 fm, how- 
ever, the full correlation function is strongly underestimated. In order to produce a real 
pion bound state one has to consider the corresponding four quark vertex and iterate this 



interaction in an RPA type approximation [16 



The zero modes are localized around the centers of the instantons with a characteristic 
size p. Evaluating the pion wave function eq. (|l|) using the propagator eq. (|i~7| ) we find 



My) = i - r7(2p) 2 + • • • (is) 



for small transverse separations y. The corresponding mean square radius is J (r 2 ) = 1.41 
fm at t = 1 fm, but does not converge as r — > oo. This means that individual instantons 
produce some interaction at short distances but the corresponding wave function has a 
long range tail and no true bound state exists. 

A similar calculation can also be done for the nucleon. Using the proton current defined 
in eq. (5) and the model for the propagator introduced above, the nucleon correlation 
function simply factorizes into the product of a free quark propagator and a scalar diquark 
correlation function. The scalar diquark correlator is very similar to the pion one, in 
particular it is also attractive and the corresponding wave function as calculated from 
the propagator eq. (fPTD is the same. The isospin structure of the proton current is such 
that the (i-quark which is used to define the Bethe Salpeter wave function is always in 
the scalar diquark. We then find that in this approximation the nucleon Bethe Salpeter 
wave function is identical to the one of the pion: iPn{v) = Vv(j/)> This fact is essentially 
a consequence of the diquark structure of the nucleon current. 

We have seen that such single-instanton approximation leads to a short range at- 
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Table 1: Rms radii extracted for the Bethe Salpeter wave functions measured at a lon- 
gitudinal separation r = 1.25 fm. The various correlation functions are described in the 

text. 





tt(P-P) 


vr(P-AV) 


p(V-V) 


p(V-T) 


\/(r 2 ) [fm] 


0.61 


0.56 


0.73 


0.70 




N (Vi - Vi) 


N(r)! - r] 2 ) 


A 




vV 2 > M 


0.63 


0.61 


0.70 





traction which affects the wave functions, although no true bound state in the pion and 
nucleon channels appears. This is even more so for the rho meson and delta channels, 
in which individual instantons do not affect the correlator at all. In order to study the 
question whether an ensemble of instantons leads to the formation of hadronic bound 
states we have done numerical simulations in the Random Instanton Model. For these 
simulations we have studied an ensemble of 256 instantons in a box 6.7 x 3.4 3 fm . In 
order to minimize finite size effects while still being able to study large separations, the 
correlation functions are measured along the long axis (the 'time' direction) of the box. 
For the measurement of wavef unctions, the available statistics limits us to the separations 
smaller than 1.5 fm. 

The results of our calculations are shown in figures 1.-4. In fig. 1. we show the measured 
pion Bethe Salpeter wave function as a function of the longitudinal separation r. The 
wave functions changes quite significantly between r = 0.5 fm and r = 1 fm, but remains 
roughly unchanged for r > 1.0 fm. We have extracted the corresponding mean square radii 
by assuming that the wave functions are exponential for transverse separations y > 1.0 fm. 



The pion rms radius measured this way is v/(r£) = 0.61 fm at r = 1.25 fm, significantly 
smaller that what we found for the dilute case considered above. 

The fact that the Bethe Salpeter wave function appears to be dominated by the pion 
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is also illustrated by a comparison of the wave functions determined from the diagonal 
and off diagonal (pseudoscalar-axialvector) correlation functions (see fig. 2. and table 
1.). The two wave functions are very similar with the off diagonal one being somewhat 
smaller. As explained in section 2 this might be due to the fact that there is still some 
contribution from the continuum left in the diagonal correlator while it is absent in the 
off diagonal one. 

In fig. 2. we also show a similar comparison for the rho meson wave function deter- 
mined from the vector and vector-tensor correlators as well as the nucleon wave function 
measured from the diagonal and the off diagonal correlator of the first and second Ioffc 
current. In both cases the agreement is excellent and we conclude that the wave functions 
are likely to be resonance dominated. This can also be seen from the fact that even in 
the rho meson channel, in which there is no first order instanton induced interaction, the 
wave function is smaller than what would be expected from the propagation of massive 
constituent quarks. 

A comparison of the wave functions obtained in different channels is shown in fig. 3. We 
observe that the pion and the proton as well as the rho meson and the delta resonance 
have very similar wave functions. The pion and the proton, however, are significantly 
smaller than the rho meson and the delta resonance. We have already argued that the 
scalar diquark content of the nucleon and the fact that instantons produce substantial 
attraction in the scalar diquark channel suggest that the nucleon and pion wave functions 
are very similar. Analogously, one can think of the delta as consisting of a quark and a 
vector diquark. In this case, however, the above argument does not really apply since the 
d quark is not necessarily sitting in the diquark and there is no interaction that would 
enhance the importance of the quark-diquark part of the correlation function. 
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In fig. 4. we compare our results for the pion and the proton with the lattice result 
reported in ||, |*|. These authors have determined the Bethe Salpeter wave functions for 
the pion, the rho meson and the nucleon in a quenched lattice simulation. In addition 
to that, they have also measured wave functions after 'cooling' their configurations. This 
procedure relaxes any given gauge field configuration to the closest classical solution 
17]. It removes essentially all gluon fluctuations except for those stabilized by topology. 



'Cooling' therefore displays the contribution of instantons to the measured wave functions. 
Comparing the different curves shown in fig. 4. we observe that the full lattice wave 
functions are more compact than the instanton results while the cooled wave functions 
are even larger. Qualitatively it is clear why the wave functions in the instanton model 
extend farther out: the instanton model misses at least two effects, the Coulomb force at 
short distances and the string potential at large distances. Both help to localize the wave 
functions, although none of them seems to be of crucial importance. 

This can be seen in more detail by considering the shape of the wave functions. At 
short distances our wave functions are quadratic in the separation, while the lattice results 
appear to be linear. This effect should be due to the perturbative Coulomb and spin-spin 
forces, since the linear behavior disappears when the wave functions are cooled. This 
can be checked further by applying only a few cooling sweeps (compared to the 25-50 
sweeps used to extract the instanton liquid): the perturbative "quantum noise" disappears 
already at this level. 

Finally, although the shapes are similar, the absolute size of the wave functions ex- 
tracted from our model and from cooled lattice configurations do not agree. A possible 
reason for this behavior is the fact that there are some ambiguities related to the redefi- 
nition of the lattice scale after cooling. Even before cooling, the lattice scale depends on 
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the particular observable that is used to fix the units [§]. 



4. Conclusions 

We have measured Bethe Salpeter wave functions for the pion, the proton, the rho 
meson and the delta resonance in the Random Instanton Liquid Model. 

In all the channels considered here there is evidence that this Model leads to the 
formation of localized wave functions. This is based on the fact that the wave functions 
become stable as the longitudinal separation is increased and the observation that the 
wave functions measured using different sources are very similar to each other. 

Furthermore, the size of these bound states are quite similar to those observed ex- 
perimentally and measured on the lattice. In fact, the model reproduces most of the 
qualitative features of the wave functions that are observed in lattice simulations. In 
particular, we find that the wave functions of the pion and the proton are very similar 
whereas the rho meson has a significantly larger size. 

The shape of the wave functions measured on the lattice are somewhat different, 
though. This is an important observation since the model does not include some of the 
physics both at very short and at very large distances, namely the perturbative and 
the confining forces. The size of the observed deviations is therefore a measure for the 
importance of these effects in the formation of hadronic bound states. After cooling the 
qualitative behavior of the lattice results agrees with what is observed in the Instanton 
Model. However, the mean square radii determined in the Instanton Model are somewhat 
smaller as compared to the cooled lattice reported in ||. 

In conclusion we have shown that the Random Instanton Liquid Model can account 
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for the qualitative behavior of hadronic Bethe Salpeter wave functions. It will therefore 
be of great interest to extend the calculations presented in this work to include the effects 
of finite temperature. 
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Figure 1: Pion Bethe Salpeter wave function measured at different separations r = 0.5. 
0.75, 1.0 and 1.25 fm. All wave functions are normalized to one at the origin. The lowest 
curve corresponds to r = 0.5 fm. The solid curves show a parametrization used to extract 
the means square radii. 



Figure 2: Bethe Salpeter wave functions for the pion, the rho meson and the nucleon 
as determined from diagonal and off diagonal correlation functions at r = 1.25 fm. The 
different correlators used are explained in the text. All wave functions are normalized to 
one at the origin. 



Figure 3: Bethe Salpeter amplitudes for the pion, the rho meson, the nucleon and the 
delta measured in the random instanton model. The wave functions are determined from 
diagonal correlation functions at r = 1.25 fm. They are normalized to one at the origin. 
The solid curves show a parametrization of the data used to extract the corresponding 
means square radii. 



Figure 4: Comparison of pion and proton Bethe Salpeter amplitudes calculated in the 
random instanton model with the corresponding results of lattice calculations reported in 
[2,3]. The lowest curve is the full lattice result, the curve in the middle shows the result 
in the instanton model and the upper curve corresponds to the lattice result after cooling. 
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1. Introduction 

Hadronic correlation functions are an important tool in understanding the complicated 
nonperturbative structure of the QCD vacuum [1]. These correlation functions are defined 
as the vacuum expectation value II(a;) =< T(jh(x)jh(0)) > of the time ordered product 
of two currents with the quantum numbers of a given hadron h. At short distances, 
asymptotic freedom predicts that the correlation function is dominated by the propagation 
of free quarks. In this case, the correlator follows a simple power law behavior, II(a;) ~ 
1/x 6 for mesons and II(a;) ~ 1/x 9 for baryons. Corrections to this behavior can be studied 
using the Operator Product Expansion. 

In general the correlation function can be expressed in terms of the contributions 
of intermediate physical states. At large spacelike separations, the lowest resonance in 
a given channel dominates and the corresponding mass and coupling constant can be 
extracted from the exponential falloff of the correlator. Roughly speaking, the coupling 
constant measures the probability for all quarks inside the hadron to be at the same 
point. In order to obtain more detailed information about the structure of the hadron 
the concept of hadronic wavefunctions was introduced in [2] . In the case of the pion this 
quantity is defined as the following gauge invariant Bethe-Salpeter amplitude 

My) = jd 4 x <0\d(x)Pe*J* +yA ^ dx ' l5 u(x + y)\TV> . (1) 

Here the two quarks are connected by a string of flux. One should note that in a relativistic 
theory the concept of a hadronic wavefunction is not uniquely defined. Other functions 
providing somewhat different information include the Bethe-Salpeter amplitudes in a fixed 
gauge or density-density correlation functions (inside a given hadron). In practice the 



Bethe-Salpeter amplitude eq. (1) is extracted from the three point correlator 

H v {x,y) = <0\T(d(x)Pe l I* +yA ^ dx ' l5 u(x + y)d(0) l5 u(0))\0>~i;(y)e- m ^ (2) 

where x has to be a large spacelike separation in order to ensure that the correlation 
function is dominated by the ground state and y is the separation of the two quarks in 
the transverse direction. 

The Bethe-Salpeter amplitude eq. (1) has been measured in a number of lattice gauge 
simulations, both at zero and at finite temperature [2, 3, 4, 5]. However, with the exception 
of one work in the context of the bag model [6] and the case of very large temperatures, 
where the wave function in the dimensionally reduced theory can be evaluated from an 
effective Schrodinger equation [7, 8], there have been no efforts to determine these wave- 
functions in a theoretical model. 

We would like to fill this gap by applying the Random Instanton Liquid Model [9, 
10, 11] to the problem of calculating hadronic wavefunctions. The model is based on 
the assumption that the QCD vacuum is dominated by a particular type of gauge field 
configurations, small size instantons. The vacuum is then characterized by the density 
n = 1 fm and the average size p = 1/3 fm of the instantons. Essentially, these numbers 
were chosen to reproduce two global properties of the QCD vacuum, the value of the 
quark and gluon condensates. Recently, we have determined a large number of hadronic 
two point correlation functions using this model [11, 12]. The model not only successfully 
describes the masses and coupling constants of the lowest resonances in most channels, 
but also reproduces the results of a recent lattice calculation of two point correlation 
functions [13]. 

In the present work we want to extend this comparison to hadronic Bethe-Salpeter 
amplitudes. First of all, this will enable us to check whether the correlators calculated 



in [11, 12] really correspond to the propagation of hadronic bound states. In addition 
to that, one can study whether the size and shape of the corresponding hadronic wave 
functions is similar to what was found in lattice simulations. In general, these are very 
non trivial questions. The Instanton Model takes into account only one very specific 
non perturbative effect among the many that determine the structure of QCD. Whether 
instanton induced effects really dominate the physics of a particular observable has to be 
checked in each individual case. 

Another interesting question concerns the differences between the wavefunctions of 
hadronic states with different quantum numbers. In a potential model approach, these 
differences are due to the gluomagnetic forces that give rise to a spin-spin interaction 
between quarks. These forces are absent in our calculations. However, there also exists 
a spin-dependent instanton induced interaction between quarks, which qualitatively de- 
scribes the mass splitting between different states. The question is whether it can also 
account for the differences in the size of the corresponding bound states. 

Finally we would like to address a few questions related to the interpretation of lattice 
results on hadronic Bethe-Salpeter amplitudes. In particular we want to study the ques- 
tions whether the observation of localized Bethe Salpeter amplitude uniquely reflects on 
the existence of a hadronic bound state and whether measurements of correlation func- 
tions at rather short distances on the order of one fermi (or even less) are sufficient to 
deduce properties of the ground state in a given channel. 



2. Wave Functions in the Random Instanton Model 



In the following we will study Bethe-Salpeter amplitudes for the pion, the rho meson, 



the nucleon and the delta resonance. For this purpose we have considered correlation 
functions involving the following currents 

j*(x) = d(x)j 5 u(x), (3) 

%{z) = d(x) 7 .u(x), (4) 

f a (x) = e abc (u a (x)C l5 d b (x))u c a (x), (5) 

*«» = z^ (2(u a (x)C^(x))ul(x) + (u a (x)Cw b (x)),r a (xj) . (6) 

Here q® denotes auor d-quark spinor with color index a and spinor index a and C is the 
charge conjugation matrix. The nucleon current given above is a linear combination of 
the standard Ioffe currents [12]. Our choice here is motivated by the fact that this current 
was used in the lattice simulations reported in [2]. We have done calculations using all 
six nucleon correlation functions introduced in [12] and will comment on the importance 
of the choice of current below. 

We consider the case of exact isospin symmetry so we have simply chosen the flavor 
content of the mesonic currents to minimize the number of contractions contributing to 
the correlation function. In the case of the nucleon we have used a proton current and 
defined the corresponding wave function by taking the two u quarks at the same point 
and measuring the correlator as a function of the separation of the d quark. Similarly, for 
the delta resonance we have taken the A + state and again measured the Bethe Salpeter 
amplitude as a function of the distance between the two u quarks and the d quark. 

The Bethe-Salpeter amplitude defined in eq. (1) involves a diagonal correlation func- 
tion of two pseudoscalar currents. The same information, however, can also be obtained 
from non diagonal correlation functions. In the case of the pion, we use the pseudoscalar- 



axialvector correlation function 

U^ aufl (x,y) = <0\T(d(x)Pe>f* +yA ^ dx ' l5 u(x + y)d(0) lfll5 u(0))\0>, (7) 

while in the case of the rho meson, we have also determined the Bethe-Salpeter amplitude 
from the off diagonal vector-tensor correlation function 

n^T^x.y) = < 0\T(d(x)Pe*J: +VA ^ dx ' 7 „u(x + y)d(0)a„„u(0))\0 > . (8) 

Similarly, for the nucleon one can use non diagonal correlation functions between the two 
different Ioffe currents [12]. The advantage of using non diagonal correlators lies in the 
fact that these correlation functions receive no contribution from free quark propagation. 
For this reason the dominance of the ground state over excited states at large but finite 
separation is expected to be somewhat better for non diagonal correlators. In any case, 
comparing the Bethe Salpeter amplitudes determined from diagonal and non diagonal 
correlators provides a nice check on whether the measured wavefunctions really do reflect 
on the groundstate in the given channel. 

In addition to the channels considered here it would also be interesting to study the 
wave functions of other mesons. Of particular interest in comparison with the rho meson 
is the ai meson. This channel, however, is difficult to analyze since the ai correlation 
function at large distances is dominated by the pion [1] and the off diagonal ai-tensor 
correlation function turns out to be very small. Furthermore, it would be nice to demon- 
strate that one does not observe localized Bethe Salpeter amplitudes in channels in which 
(at least inside the model) no bound states are expected to form. In the Instanton Model, 
one such channel would be the scalar isovector (£-meson) channel. Again, the analysis 
turns out to be very difficult since the fact that the interaction is repulsive implies that the 
measured correlation function at large distances is very small so that the wave function 
cannot be reliably determined. 



The instanton ensemble in QCD can be generally described by a partition function of 
the type 

Z = / U l [dQ, l exp(-S l )] exp(-S int )Ii.f =1)Nf det(ii) + im f ) (9) 

where [dQi] denotes the integration over the collective coordinates (orientation, size, posi- 
tion) associated with the instantons, Si = 8tt 2 /g 2 (p) is the action of an individual instan- 
ton and Si n t is the instanton-instanton interaction. Since the evaluation of this partition 
function still constitutes a formidable problem we consider the even simpler Random In- 
stanton Model. In this model we assume that except for the size, which we keep fixed, the 
distribution of the collective coordinates is completely random. In particular, adopting 
this assumption we neglect the effects of the fermion determinant det(iD + imf), which 
is similar to "quenched" lattice simulations we are going to compare with. 

The Bethe-Salpeter amplitudes are then determined from the quark propagator in a 
given configuration. The pion correlation function eq. (2), for example, is given by 

IL^(x,y) = <TiS ab (x + y)(Pe l f* +yA{x ' ]dx ') bc S ca (-x)> (10) 

where the trace is performed with respect to the Dirac indices. The calculation of the 
propagator is explained in detail in [11]. Here we only outline the general strategy. The 
quark propagator in the field of a single instanton can be written as the sum of a zero 
mode contribution and the propagator due to non zero modes. Both pieces are known 
analytically. In the many instanton configuration we treat the zero modes exactly by 
numerically inverting the Dirac operator in the space spanned by the zero modes of the 
individual instantons. The non zero modes, on the other hand, are taken into account 
only approximately, by including the first term in a multiple scattering expansion. The 



resulting quark propagator is given by 

s ( x ,y) = J2 < x \^i >< ^i\~ — f; — -^\4>j >< 4>j\y > +S N zM(x,y), (n) 

where (j>i(x) is the zero mode associated with the instanton / and 

Snzm(x,v) = S (x,y) + Ys( S NZM( x ,y) ~ S (x,y)) (12) 

is the propagator due to non zero modes. Here, So denotes the free propagator and S T NZM 
is the non zero mode propagator in the field of an individual instanton. The non zero 
mode propagator also has to be corrected for the effects of current quark masses [11]. In 
this paper we take the light quark mass to be m Uj( j = 10 MeV and the corresponding 
effects are small. The Schwinger factor P exp(i / A^dx^) is calculated numerically for any 
given gauge field configuration by expanding the pathordered exponential as an infinite 
product 

Pexp(i / Af.dXf,) = 11,(1 + iA^x^dx^) . (13) 

In practice this expression is approximated by a finite product where the stepsize dxi is 
determined by the magnitude of the local gauge potential. 

In general, the inclusion of the Schwinger factor eq. (13) is expected to give an impor- 
tant contribution to the measured wave functions, since it corresponds to an additional 
string type potential. The instanton liquid, however, does not produce confinement and 
no string potential is expected to appear. At large separations, the Schwinger factor 
decreases exponentially Pexp(i / r Adx) ~ exp(— mr) with m being the mass renormal- 
ization of a heavy quark in the instanton background [12]. Using the 'sum ansatz' for 
the gauge field of the multi-instanton configuration we find to = 65 MeV, in agreement 
with estimates made in [14]. Employing the more realistic ratio ansatz an even smaller 



screening mass to = 20 MeV is found. In any case, the insertion of the path ordered 
exponential will only lead to minor effects on the wave functions we are going to discuss. 

3. Results and comparison with Lattice Gauge Simulations 

Before we come to the results obtained, let us discuss the Bethe Salpeter wave functions 
obtained by using simplified models for the propagator. Since in all numerical simulations 
to date the separation between the sources is rather modest (at most on the order of 1 
fm) it is important to understand what the size of non asymptotic effects in the measured 
wave functions are. 

In the simplest case we consider a free fermion propagator 

S 'W = i 2 ^ < 14 > 

and take the path ordered exponential to be unity. Here all gamma matrices and coor- 
dinates are taken to be euclidean. Inserting this propagator in the pion correlator eq.(2) 
we get 1 

Kf 3. 2 _ ]P_ N c 2 

At distances r = \Jx 2 on the order of one fermi, this corresponds to a 'wave function' 



with a mean square radius 2 J (r 2 ) = 0.89 fm, quite different from the asymptotic result. 
This effect becomes even more pronounced if one considers the propagation of a massive 
constituent quark 

T * X 

S(x)=i D'(m q ,T) + m q D(m q ,T) (16) 



1 Since we are discussing short range effects, we are using a symmetric configuration in which the 
quarks are located at transverse separations —y/2, +y/2. For asymptotic wave functions, this difference 
should not matter. 

2 The mean square radius is defined by (r 2 ) = J d 3 rr 2 \tp(r)\ 2 where the Bethe Salpeter wave function 
is normalized to one. 



where D(m,T) = m/(4:Tr 2 T)Ki(mT) is the coordinate space propagator of a massive scalar 
particle. Taking the quark mass to be 300 MeV, we find a wave function with a size of 



(r 2 ) = 0.79 fm at a distance of r = 1.0 fm. Even at a distance of 1.5 fm, the wave 



function is still fairly well localized with a radius of J (r 2 ) = 1.1 fm. Note however, that 
so far the appearance of a localized wave function is a purely geometric effect and has 
nothing to do with the existence of bound states. 

One should note that the fact that the Bethe Salpeter wavefunctions have not reached 
their asymptotic shape can not be detected from an analysis of the spacelike screening 
mass. The screening masses in this simple model are m scr = 2m q for mesons and m scr = 
3m q for baryons. There are only small power corrections to this behavior that are due 
the fact that the quarks can have some transverse momentum. We would also like to 
emphasize that considering the effects of finite quark masses on the Bethe Salpeter wave 
functions is very important at finite temperature where quarks develop thermal screening 
masses m scr = ttT. Around T c this mass is on the order of 500 MeV. If the wave functions 
are determined at distances x ~ 1/T, which is typical for current lattice simulations, 
masses will have an important effect on the measured Bethe Salpeter wave functions. 

After these remarks concerning non asymptotic effects we would like to study how 
instantons lead to the formation of bound states that manifest themselves in localized 
asymptotic Bethe Salpeter amplitudes. A simple one-instanton approximation that has 
been used [15] to analyze the qualitative effects of instantons on the short- distance be- 
haviour of correlation functions is based on the following expression for the propagator in 
a dilute instanton gas 

S(x) = b (x) + n dz — 7—^ . (17) 

J blYh 

Here n is the density of instantons, (j) 0jZ (x) is a zero mode localized at z and m* = 



2Tvp c Jn /Q is the effective quark mass in the instanton gas. Recall, that such propagator 
eq. (17) leads to an attractive pion correlation function and explains the failure of the 
operator product expansion in this channel. At distances greater than r = 0.5 fm, how- 
ever, the full correlation function is strongly underestimated. In order to produce a real 
pion bound state one has to consider the corresponding four quark vertex and iterate this 
interaction in an RPA type approximation [16]. 

The zero modes are localized around the centers of the instantons with a characteristic 
size p. Evaluating the pion wave function eq. (1) using the propagator eq. (17) we find 

Mv) = i - y 2 /( 2 p) 2 + • • • (is) 



for small transverse separations y. The corresponding mean square radius is J (r 2 ) = 1.41 
fm at t = 1 fm, but does not converge as r — > oo. This means that individual instantons 
produce some interaction at short distances but the corresponding wave function has a 
long range tail and no true bound state exists. 

A similar calculation can also be done for the nucleon. Using the proton current defined 
in eq. (5) and the model for the propagator introduced above, the nucleon correlation 
function simply factorizes into the product of a free quark propagator and a scalar diquark 
correlation function. The scalar diquark correlator is very similar to the pion one, in 
particular it is also attractive and the corresponding wave function as calculated from 
the propagator eq. (17) is the same. The isospin structure of the proton current is such 
that the d-quark which is used to define the Bethe Salpeter wave function is always in 
the scalar diquark. We then find that in this approximation the nucleon Bethe Salpeter 
wave function is identical to the one of the pion: ijj^{y) = ^(y). This fact is essentially 
a consequence of the diquark structure of the nucleon current. 

We have seen that such single-instanton approximation leads to a short range at- 
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Table 1: Rms radii extracted for the Bethe Salpeter wave functions measured at a lon- 
gitudinal separation r = 1.25 fm. The various correlation functions are described in the 
text. 





tt(P-P) 


tt(P-AV) 


KV-V) 


P (Y-T) 


\/(r 2 ) M 


0.61 


0.56 


0.73 


0.70 




N( Vl - Vl ) 


N( Vl - 7/ 2 ) 


A 




\/(r 2 ) M 


0.63 


0.61 


0.70 





traction which affects the wave functions, although no true bound state in the pion and 
nucleon channels appears. This is even more so for the rho meson and delta channels, 
in which individual instantons do not affect the correlator at all. In order to study the 
question whether an ensemble of instantons leads to the formation of hadronic bound 
states we have done numerical simulations in the Random Instanton Model. For these 
simulations we have studied an ensemble of 256 instantons in a box 6.7 X 3.4 3 fm . In 
order to minimize finite size effects while still being able to study large separations, the 
correlation functions are measured along the long axis (the 'time' direction) of the box. 
For the measurement of wavefunctions, the available statistics limits us to the separations 
smaller than 1.5 fm. 

The results of our calculations are shown in figures 1.-4. In fig. 1. we show the measured 
pion Bethe Salpeter wave function as a function of the longitudinal separation r. The 
wave functions changes quite significantly between r = 0.5 fm and r = 1 fm, but remains 
roughly unchanged for r > 1.0 fm. We have extracted the corresponding mean square radii 
by assuming that the wave functions are exponential for transverse separations y > 1.0 fm. 



The pion rms radius measured this way is J (r^) = 0.61 fm at r = 1.25 fm, significantly 
smaller that what we found for the dilute case considered above. 

The fact that the Bethe Salpeter wave function appears to be dominated by the pion 
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is also illustrated by a comparison of the wave functions determined from the diagonal 
and off diagonal (pseudoscalar-axialvector) correlation functions (see fig. 2. and table 
1.). The two wave functions are very similar with the off diagonal one being somewhat 
smaller. As explained in section 2 this might be due to the fact that there is still some 
contribution from the continuum left in the diagonal correlator while it is absent in the 
off diagonal one. 

In fig. 2. we also show a similar comparison for the rho meson wave function deter- 
mined from the vector and vector-tensor correlators as well as the nucleon wave function 
measured from the diagonal and the off diagonal correlator of the first and second Ioffe 
current. In both cases the agreement is excellent and we conclude that the wave functions 
are likely to be resonance dominated. This can also be seen from the fact that even in 
the rho meson channel, in which there is no first order instanton induced interaction, the 
wave function is smaller than what would be expected from the propagation of massive 
constituent quarks. 

A comparison of the wave functions obtained in different channels is shown in fig. 3. We 
observe that the pion and the proton as well as the rho meson and the delta resonance 
have very similar wave functions. The pion and the proton, however, are significantly 
smaller than the rho meson and the delta resonance. We have already argued that the 
scalar diquark content of the nucleon and the fact that instantons produce substantial 
attraction in the scalar diquark channel suggest that the nucleon and pion wave functions 
are very similar. Analogously, one can think of the delta as consisting of a quark and a 
vector diquark. In this case, however, the above argument does not really apply since the 
d quark is not necessarily sitting in the diquark and there is no interaction that would 
enhance the importance of the quark-diquark part of the correlation function. 
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In fig. 4. we compare our results for the pion and the proton with the lattice result 
reported in [2, 3]. These authors have determined the Bethe Salpeter wave functions for 
the pion, the rho meson and the nucleon in a quenched lattice simulation. In addition 
to that, they have also measured wave functions after 'cooling' their configurations. This 
procedure relaxes any given gauge field configuration to the closest classical solution 
[17]. It removes essentially all gluon fluctuations except for those stabilized by topology. 
'Cooling' therefore displays the contribution of instantons to the measured wave functions. 
Comparing the different curves shown in fig. 4. we observe that the full lattice wave 
functions are more compact than the instanton results while the cooled wave functions 
are even larger. Qualitatively it is clear why the wave functions in the instanton model 
extend farther out: the instanton model misses at least two effects, the Coulomb force at 
short distances and the string potential at large distances. Both help to localize the wave 
functions, although none of them seems to be of crucial importance. 

This can be seen in more detail by considering the shape of the wave functions. At 
short distances our wave functions are quadratic in the separation, while the lattice results 
appear to be linear. This effect should be due to the perturbative Coulomb and spin-spin 
forces, since the linear behavior disappears when the wave functions are cooled. This 
can be checked further by applying only a few cooling sweeps (compared to the 25-50 
sweeps used to extract the instanton liquid): the perturbative "quantum noise" disappears 
already at this level. 

Finally, although the shapes are similar, the absolute size of the wave functions ex- 
tracted from our model and from cooled lattice configurations do not agree. A possible 
reason for this behavior is the fact that there are some ambiguities related to the redefi- 
nition of the lattice scale after cooling. Even before cooling, the lattice scale depends on 
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the particular observable that is used to fix the units [2] 



4. Conclusions 

We have measured Bethe Salpeter wave functions for the pion, the proton, the rho 
meson and the delta resonance in the Random Instanton Liquid Model. 

In all the channels considered here there is evidence that this Model leads to the 
formation of localized wave functions. This is based on the fact that the wave functions 
become stable as the longitudinal separation is increased and the observation that the 
wave functions measured using different sources are very similar to each other. 

Furthermore, the size of these bound states are quite similar to those observed ex- 
perimentally and measured on the lattice. In fact, the model reproduces most of the 
qualitative features of the wave functions that are observed in lattice simulations. In 
particular, we find that the wave functions of the pion and the proton are very similar 
whereas the rho meson has a significantly larger size. 

The shape of the wave functions measured on the lattice are somewhat different, 
though. This is an important observation since the model does not include some of the 
physics both at very short and at very large distances, namely the perturbative and 
the confining forces. The size of the observed deviations is therefore a measure for the 
importance of these effects in the formation of hadronic bound states. After cooling the 
qualitative behavior of the lattice results agrees with what is observed in the Instanton 
Model. However, the mean square radii determined in the Instanton Model are somewhat 
smaller as compared to the cooled lattice reported in [3]. 

In conclusion we have shown that the Random Instanton Liquid Model can account 
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for the qualitative behavior of hadronic Bethe Salpeter wave functions. It will therefore 
be of great interest to extend the calculations presented in this work to include the effects 
of finite temperature. 
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Figure 1: Pion Bethe Salpeter wave function measured at different separations r = 0.5, 
0.75, 1.0 and 1.25 fm. All wave functions are normalized to one at the origin. The lowest 
curve corresponds to r = 0.5 fm. The solid curves show a parametrization used to extract 
the means square radii. 
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Figure 2: Bethe Salpeter wave functions for the pion, the rho meson and the nucleon 
as determined from diagonal and off diagonal correlation functions at r = 1.25 fm. The 
different correlators used are explained in the text. All wave functions are normalized to 
one at the origin. 
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Figure 3: Bethe Salpeter amplitudes for the pion, the rho meson, the nucleon and the 
delta measured in the random instanton model. The wave functions are determined from 
diagonal correlation functions at r = 1.25 fm. They are normalized to one at the origin. 
The solid curves show a parametrization of the data used to extract the corresponding 
means square radii. 
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Figure 4: Comparison of pion and proton Bethe Salpeter amplitudes calculated in the 
random instanton model with the corresponding results of lattice calculations reported in 
[2,3]. The lowest curve is the full lattice result, the curve in the middle shows the result 
in the instanton model and the upper curve corresponds to the lattice result after cooling. 
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